FISIST/ll-99/CFIF 



|hep-ph/9907399| 



Resonance Spin Flavour Precession and Solar 

Neutrinos 



Joao Pulido Q and E. Kh. Akhmedov f\ 

Centro de Fisica das Interacgdes Fundamentals ( CFIF) 
Departamento de Fisica, Instituto Superior Tecnico 
Av. Rovisco Pais, P-1049-001 Lisboa, Portugal 



Abstract 

We examine the prospects for the resonance spin flavour precession as a solution to 
the solar neutrino problem. We study seven different realistic solar magnetic field pro- 
files and, by numerically integrating the evolution equations, perform a fit of the event 
rates for the three types of solar neutrino experiments (Ga, CI and SuperKamiokande) 
and a fit of the energy spectrum of the recoil electrons in SuperKamiokande. A x 2 
analysis shows that the quality of the rate fits is excellent for two of the field profiles 
and good for all others with x 2 /d.o.f. always well below unity. Regarding the fits 
for the energy spectrum, their quality is better than that for the small mixing angle 
MSW solution of the solar neutrino problem, at the same level as that for the large 
mixing angle MSW solution but worse than that for the vacuum oscillations one. The 
experimental data on the spectrum are however largely uncertain especially in the 
high energy sector, so that it is too early yet to draw any clear conclusions on the 
likeliest type of particle physics solution to the solar neutrino problem. 
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1 Introduction 



More than thirty years after its first recognition [lj], the solar neutrino problem has become 
well established on the grounds of theoretical [p|-|J10| versus experimental developments 
11|-[14|. The neutrino oscillation solutions to the anomaly, both matter |1| and vacuum 
oscillations [16], have received a great deal of attention [17]-[24] and became the most 
popular ones. On the other hand, the solution based on the resonance spin flavour precession 
(RSFP) mechanism, proposed in 1988 [2J|, involves the simultaneous flip of both neutrino 
chirality and flavour [26). It has been much less thoroughly investigated, in particular 
concerning its parameter range predictions |27|] and the quality of the fits involved. In 
fact, it is only recently that the first detailed data analysis based on the minimum y 2 
approach has been performed in the framework of this mechanism [28|. One of the reasons 



for this comparative lack of popularity may be the large order of magnitude of the neutrino 
magnetic moment required for a significant conversion and our relatively limited knowledge 
of the solar magnetic field. It should however be noted that models for a large neutrino 
magnetic moment not conflicting with the smallness of neutrino mass already exist [29], and 
there are solar model indications favouring a field as large as 3 x 10 5 G around the bottom of 

As shall be seen, the required neutrino magnetic moment is, on the 



the convective zone 30 



other hand, not inconsistent with astrophysical limits |31|, especially taking into account 
the inherent uncertainties of the analyses. 

One of the main motivations for investigating the RSFP mechanism is that neutrinos 
provide a unique probe for the interior solar magnetic field, if they have a sizeable magnetic 
moment. In fact, previous analyses [27], [32]-[34] favour solar magnetic fields rising along a 
relatively short distance (a factor of at least 6-7 over a 7-9% fraction of the solar radius at 
the most). This rise may be very sharp and even discontinuous profiles are allowed by the 
data. This appears as a natural consequence of the strong suppression for the intermediate 
energy neutrinos ( 7 Be and CNO fluxes) together with the almost no suppression for the low 
energy pp ones. Solar physics arguments suggest that such a sharp rise, if it exists, must 
lie around the upper layers of the radiative zone and the bottom of the convective zone. 
Moreover the moderate suppression of the energetic 8 B neutrinos favours a gradual decrease 
of the field along the convective zone, with a larger slope at a greater depth. 

In this paper we perform an investigation of the RSFP scenario for seven different 
magnetic field profiles all obeying the general features described above. We start with the 
numerical integration of the evolution equations for two neutrino flavours, from which we 
obtain the survival probability as a function of energy. This is used to obtain the event rate 
in each of the three types of experiments (Gallium, Chlorine and SuperKamiokande). The 
rates thus obtained are then confronted with the experimental ones jn}|-||14||. We let the 
parameters Am^ (mass square difference between neutrino flavours) and Bq (value of the 
field at the peak) run over their "plausible" ranges, namely 4 x 10~ 9 eV 2 < Am^ ^ 20 x 10 -9 
eV 2 and 3 x 10 4 G < B < 3 x 10 5 G. These ranges are "plausible" in the sense that they 
are dictated by consistency with the RSFP scenario. The first is required by the location of 
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neutrino resonances as a function of energy, so that the most strongly suppressed neutrinos 
have their resonance around the strongest magnetic field, and the second range by the solar 
physics requirement of keeping Bq not greater than 3 x 10 5 G at the bottom of the convective 
zone but large enough to produce significant conversion. For each magnetic field profile we 
select those intervals in which the theoretical ratios approach their experimental values, and 
perform a x 2 analysis in order to determine the local minima and the best fit. We thus get 
impressively small values of x 2 /d.o.f., in the range (3 — 5) x 10~ 2 at the best fit in three of 
the profiles, two of which are remarkably stable against variations of the parameters Am^ 
and Bq. We also compare the theoretical prediction for the recoil electron energy spectrum 
in SuperKamiokande in each best fit case with the experimental result. 

The paper is organized as follows: in section 2 we present the magnetic field profiles 
to be investigated and the survival probabilities obtained from them upon integration of 
the evolution equation. In section 3 we calculate the ratios of the event rates predicted in 
the RSFP scenario to the standard solar model (SSM) event rate in each experiment and 
analyse their dependence on B and Am^. Next we describe our x 2 analysis for the rates 
and identify the local and absolute minima in each field profile in terms of these parameters. 
The x 2 analysis for the electron recoil spectra in SuperKamiokande ends this section. Finally 
in section 4 we present our summary and overview. 



2 Magnetic Field Profiles and Survival Probabilities 



In this section we will present the magnetic field profiles to be used throughout, all satisfying 
the general features explained in the introduction, with a sudden rise around the bottom of 
the convective zone, at approximately 0.71 of the solar radius, and a smoother decrease up 
to the surface. Here the field intensity should not exceed a few hundred Gauss. All profiles 
analysed are displayed in figs. 1 and 2. The first two are simple triangle profiles, with the 
field intensity rising linearly from zero at fraction x = xr of the solar radius, reaching a 
peak at x = xc and decreasing linearly to zero at the solar surface [B3 



B = B, 



B = , x < x R (1) 

X — Xr 

o , x R <x <x c (2) 

X C - Xr 

B = B Q \l- ^—^1 , x c < x < 1 (3) 
L 1 — x c i 

with units in Gauss. We study this class of profiles in two cases shown in fig. 1: profile 1, 
xr = 0.70, xc = 0.85 (equilateral triangle) and profile 2, xr = 0.65, xc = 0.80 Q. 

3 A slight variant of this second profile with a vanishing field for x > 0.95 was used in ref. |28| . 
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Next we consider a slightly larger upgoing slope and a quadratic like decrease up to the 
surface. This is profile 3 with xr = 0.65, xc = 0.75 (fig. 1) 



B = B, 



B = , 

X — Xr 



0" 



B = Bn 



X C - Xr 

'x-0.7\ 2 ' 



X < Xr 
Xr<X<X C 



0.3 



, Xc < x < 1. 



(4) 
(5) 

(6) 



The two cases with an infinite slope at x = xr are considered next. The first is profile 4 
with the same quadratic like decrease from x > xr up to the surface [[32] (fig. 1) 



B = Bn 



B = , X < Xr 

' .x-0.7^ 21 



0.3 



X > Xr, 



(7) 
(8) 



with xr = 0.71 and the second is profile 5 with a reversed shape slope: a strong decrease 
at first, which then becomes quite moderate on approaching the surface [B3, EM (fig. 2) 



B = , X < Xr (9) 

B = r , n f ° r,*>a* (10) 

cosh 30 (x — Xij) 

with 1^ = 0.71. Profile 6 is a modification of the previous one with a large finite positive 
slope across the bottom of the convective zone |27j (fig. 2) 



B = B X 
B = 



£ = 2.16 x 10 3 

x -0.75\ 2 
0.04 J 

B 



x < 0.7105 
0.7105 < x < 0.7483 

, 0.7483 < x < 1 



(11) 
(12) 

(13) 



cosh 30(x- 0.7483) 

with Bo = 0.998-Bi. The last case considered is profile 7, similar to the previous one but 
with a linear decrease towards the surface |27|] (fig. 2) 



B = 2.16 x 10 J 



B = B () 



1 



x-0.75\ 2 
0.04 ) 



x < 0.7105 
0.7105 < x < 0.7483 



(14) 
(15) 



B = 1.149450(1 - 3.4412(a; - 0.71)] , 0.7483 < x < 1. (16) 

The reasons to select these classes of solar field profiles are twofold: first to investigate 
previously proposed cases in the light of the most recent experimental data and second to 
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explore those profiles which provide the best and most stable fits against variations of Am|i 
and Bq. While, as will be seen, profile 3 provides the best fit of all, this is only marginally 
better than the fits corresponding to profiles 5 and 6 which are on the other hand much 
more stable, clearly favouring a specific class of field shape. 

We consider the neutrino propagation equation through solar matter which for the two 



generation Majorana case and zero vacuum mixing angle reads [35, 36 



l dt{Z L R ) [hIb (1T) 

Here V e and are the matter induced potentials for electron neutrinos and muon antineu- 
trinos respectively, and fi u is the neutrino transition magnetic moment. 

For definiteness we consider v eL — > v iiR transitions but all our results apply to v eL — > 
v TR transitions as well. We then proceed with the integration of this system of equations 
using the solar density and neutrino fluxes |57|, thus obtaining the survival and conversion 



probabilities of left handed electron neutrinos. In our numerical calculations we use \i v = 
10 _11 //b. Since the neutrino magnetic moment enters into the evolution equation only in 
the combination fi u B, the results would also apply to any other value of fi u provided that 
the solar magnetic field strength is rescaled accordingly. For each magnetic field profile the 
survival probability is given in figs. 3 and 4 for the values of Am^ and B corresponding 
to the respective best fit. We note the almost total survival of low energy neutrinos, the 
strong suppression of intermediate energy ones and the moderate suppression of the most 
energetic 8 B ones. The high energy limit is, in the RSFP case, close to 1/2, as confirmed 
by these figures, a different situation from the small mixing angle MSW where this limit is 
close to unity. This fact leads in general to better fits of the rates in the case of the RSFP 
mechanism than in the MSW one. This is because in the case of the MSW mechanism 
one has to choose very carefully the value of Am^ in order to achieve a factor of ~ 1/2 
reduction of the high energy portion of 8 B neutrinos, whereas this comes out automatically 
in the case of the RSFP mechanism. 

The next step is to insert the survival probability into the expression for the event rate 
in each experiment and take the ratio to the SSM event rate. The event rate of the gallium 
experiments in the framework of the RSFP scenario is calculated as 

RGa = ^2 R-Ga,i p Ga,t (18) 

where the sum extends over the relevant neutrino fluxes, (i = pp, pep, 7 Be, 8 B, 13 N, 15 0). 
The quantities Raa,i denote the partial event rate for flux i, as predicted by the RSFP 
model, where, except for % = 7 Be, pep we have 

R Ga ,i = / a Ga (E)P(E)f t (E)dE. (19) 



m i n 
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Here fi(E) represents the i-ih flux from ref. |37|], and the integral is taken from the experi 



mental threshold E min = E t h = 0.236 MeV for the pp, 15 0, 13 N and 8 B cases. For i = 7 Be, 
pep, the quantities Rc a ,i are just reduced to the survival probability at the corresponding 
energy. The quantities Rc a ,i are on the other hand the corresponding SSM predictions for 
the partial rates |§ obtained from (19) with P = 1 (see also table I). The RSFP predicted 
ratio of the rates for the Gallium experiment is finally 

rca = |^ (20) 
where Rc a is the SSM total rate for the gallium experiment, Rc a = J2i RGa,i- 



Source 


Flux 




Ga 


CI 




(10 10 cm- 2 s 




(SNU) 


(SNU) 


PP 


5.94 




69.6 


0.0 


pep 


1.39 x 10~ 


-2 


2.8 


0.2 


hep 


2.1 x 10" 


7 


0.0 


0.0 


7 Be 


4.8 x 10" 


1 


34.4 


1.15 


8 B 


5.15 x 10" 


-4 


12.4 


5.9 


13 N 


6.05 x 10" 


-2 


3.7 


0.1 


15Q 


5.32 x 10- 


-2 


6.0 


0.4 


17 F 


6.33 x 10- 


-4 


0.1 


0.0 


Total 






129± 8 


7 7± L2 

' • ' ^l.O 



Table I - Fluxes and partial rates (Rji coefficients in the main text) as predicted by the 
solar standard model 0. 

In order to obtain rci this procedure is repeated for the chlorine experiment with the 
obvious omission of pp neutrinos and the truncation of the other continuous fluxes below 
E t h = 0.814 MeV, the energy threshold in the CI experiment. The SSM rates Rci,i are also 
listed in table I. Regarding SuperKamiokande only one neutrino flux is present, namely 8 B 
0, and we have for the RSFP predicted ratio 

R _ it:: me) it^pj^M + 1 1 - p ( E )^ dTdE m) 
x£r me) jgr ^dTdE 



Here T is the electron recoil energy, T min = E e±h — m e with E e±h = 5.5 MeV [14| 



4 The contribution of a small flux of hep neutrinos is completely negligible in the total rates but may be 
of some importance for the high energy part of the recoil electron spectrum in SuperKamiokande. We shall 
discuss this point in more detail in section 4. 
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and the integration limits for the neutrino energy are respectively, 



T • -I- \/T 2 ■ 4- 2m T ■ 

77, ± mm ~ V ± mm \ ^' ,b e ± mm /'OQ^ 
&min — ^ V ) 

and E max = 15 MeV. The weak interaction differential cross section for v e e scattering is 
given by 

^ Gl?2 \{9v + 9a) 2 + (9v - 9a) 2 (l-l) 2 - (9v 2 - 9A 2 )^} (24) 



dTdE 2vrm e L ^ v ' » AJ ' ^ ^ V EJ ^ v ^ ' E 2 



with gy = \ + 2sm 2 6*vi/, 9a—\- The cross section for z/^e and i/ r e scattering is on the other 
hand 

with g v = — ~ + 2sin 2 9w, gA = — §■ The contribution of the neutrino electromagnetic cross 
section to eq. (21) is negligible, typically 3 orders of magnitude lower than the weak cross 
section and will be omitted. We do not include the electron energy resolution function of 
the detector as its effect on the total rate is negligible; we will, however, include it in the 
calculation of the recoil electron spectrum. We calculate the ratios r Ga ,r C i, Rsk m the 
parameter ranges Am|i = (4 — 20) x 1CT 9 eV 2 and B = (3 — 30) x 10 4 G for all seven 
magnetic field profiles presented above and identify in each case those intervals in which all 
three ratios simultaneously approach their experimental values (table II). We then perform 
for each interval a x 2 analysis in order to evaluate the quality of the fit. 



Experiment 


Data 


Theory 


Data/Theory 


Reference 


Homestake 


2.56 ±0.16 ±0.15 


7 7+ 1 -' 2 


0.33 ±0.03 




11 




Gallex 


76.4 ± 6.3±H 


129±§ 


0.59 ±0.06 


|13j 




SAGE 


69.9± 8 7 ;?±i;?' 


129±^ 


0.54 ±0.06 




VA 




SuperKamiokande 


2.44 ± 0.05±g;g§ 


5.i5±j;? 


0.474 ±0.020 


\U\ 



Table II - Data from the four solar neutrino experiments. Units are SNU for the first 
three experiments and in cm~ 2 s~ l for SuperKamiokande. For the Gallium experiments we 
have used the combined result 72.3 ± 5.6 SNU. 



3 Analysis of the data 
3.1 Description of x 2 analysis 

We follow in this subsection the procedure outlined in ref. |38| adapted to the RSFP 



mechanism and the present experimental situation |IH-|I4]|. We start with the definition of 
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the x 2 function 



3 - -i 



X 2 = E (Rn - Rn XP ) [* 2 (tot) _ (R h - R n ^). (26) 
iu'2=i 



JU2 



In this expression the indices j\iH run over the three types of experiments: Ga, CI and 
SuperKamiokande respectively. The quantities R denote the theoretical (RSFP) event rates 
defined in eqs. (18), (19) for Ga and CI experiments, and in (20) for the SuperKamiokande 
one and the quantities R exp denote the experimental ones given in table II. They are ex- 
pressed in SNU for Ga, CI and in terms of the ratio data/SSM for SuperKamiokande. For 
the Gallium case we use the combined value from SAGE and Gallex, R^ = 72.3 ±5.6 SNU. 
The matrix c 2 j 1 j 2 (tot) is the total error matrix given by 

a 2 jlj2 (tot) = a 2 jlj2 (exp) + a 2 jlj2 (th) . (27) 

The experimental error matrix is given in terms of o~j exp , the experimental errors shown in 
table II, 

a 2 nn (exp) = 5 nn a n exp v» exp (28) 

so that 

/ 31.36 \ 

cr 2 jlj2 (exp) = 0.0481 . (29) 
V 0.0004 / 

The theoretical error matrix in eq. (27) is the sum of the cross section contribution con- 
taining the uncertainties in the cross sections with the astrophysical contribution involving 
the uncertainties in the astrophysical parameters 

^hh ( th ) = ^jih M + ° 2 hh M • ( 30 ) 

The uncertainties in the detector cross sections, expressed in the form AlnCV,', entering the 
cross section error matrix, are assumed uncorrelated and were taken from ref. 0]. They 
give rise to the following diagonal matrix 

4i» = ^E#V(AlnC^) 2 . (31) 

i=\ 

where the index i runs over all relevant fluxes, i = pp, pep, 7 Be, 8 B, 13 N, 15 0. Model 
dependence appears in this matrix through the R'^s. 

The astrophysical error matrix is finally 

6 9 

° 2 hhM= E R hh R i2j2 J2a hk a i2k (A\nX k ) 2 (32) 

81,22=1 k=l 

In this expression index k extends over the nine astrophysical parameters k = Sn, S 33 , 634, 
-Si, 14, 5*17, Lum, Z/X, Age, Opac, the a matrix denotes the logarithmic derivatives of the six 
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fluxes considered with respect to these parameters, 
notation is self clear. The (6 x 9) matrix 

a ik A In X k 



aik = dhi fi/dhiXk and the rest of the 



(33) 



which appears in equation (32) was calculated using the data from f37[. We now have all 
the necessary elements to calculate y 2 for each magnetic field profile. 



3.2 Fits of the rates and their quality 

In this subsection we will investigate the local minima of x 2 for all studied solar field profiles, 
identifying in each profile the lowest minimum in terms of B (peak field value) and Am^. 
As discussed in the introduction we analyze the intervals Bq = (3 — 30) x 10 4 G and 
Am^i = (4 — 20) x 10~ 9 eV 2 , the first interval dictated by the requirement of significant 
neutrino conversion (lower bound) and solar physics arguments [[3(J (upper bound) and 
the second interval dictated by the resonance locations of the several neutrino fluxes. In 
fitting the rates we have three types of experiments (Ga, CI, SuperKamiokande) and two 
parameters: Am^, [i v Bq for a given magnetic field profile. Hence the number of degrees 
of freedom is one. We search for values of \ 2 smaller than or of order unity and discuss 
the profiles in the same order as in section 2. Starting therefore with profile 1 (equilateral 
triangle), eqs. (l)-(3), we observe two local minima (fig. 5), one at B ~ 1.68 x 10 5 G and 
the other at B ~ 2.71 x 10 5 G, both with Am^ ~8x 10~ 9 eV 2 . The first is an absolute 
minimum with x 2 m in = 0.085 and for the second, x 2 m in = 0.137. While such values of 
X 2 (especially the first one) are remarkably low, the stability of both fits against small 
parameter variations is rather poor: a change of 3 kG in the value of Bq implies a change 
in x 2 by a factor of order 30. 

The second profile is the triangle field, eqs (l)-(3), in which we identified three local 
minima (fig. 6): the best fit {x 2 m in — 0.10) corresponds to B ~ 1.23 x 10 5 G, Am^ ~ 
1.20 x 10~ 8 eV 2 , the second best (x 2 min = 0.284) to B ~ 1.98 x 10 5 G, Am 2 21 ~ 1.20 x 10~ 8 
eV 2 and the third (x 2 min = 0.553) to B ~ 2.85 x 10 5 G, Am 2 21 ~ 1.25 x 10" 8 eV 2 . These 
minima are displayed in fig. 6 and it is seen that their stability is better than the previous 
cases: a change of 3 kG implies a change in x 2 by a factor of 3-4. 

Next we examine profile 3, eqs. (4)- (6): there are two local minima of x 2 both corre- 
sponding to Am^ ~ 1.2 x 10 -8 eV 2 (fig. 7), the first one at B = 9.54 x 10 4 G with an 
impressively low x 2 (Xmin = 0.036) and the second one at B ~ 1.72 x 10 5 G (x 2 m in = 0.11). 
We note the extreme instability of the first, where the same change of 3 kG in Bq leads to a 
change in x 2 by a factor of 50, while for the second the corresponding factor is of order 30. 

Profile 4, eqs. (7), (8), shows two local minima of x 2 (fig- 8) : X 2 m in = 0-68 for 
Am 2 ! ~ 1.2 x 10~ 8 eV 2 , B = 9.0 x 10 4 G and x 2 min = 0-59 for Am 2 21 ~ 1.6 x 10~ 8 eV 2 , 
B = 1.7 x 10 5 G. These fits are poorer than the previous ones and in addition they are also 
quite unstable: a factor of at least 10 increase in x 2 f° r a 3 kG change in the field at the 
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peak. Of all cases examined, such a profile appears to be the least favoured by the data. If 
in this case we choose a power greater than 2 in eq. (8), the results become even worse as 
the SuperKamiokande and Chlorine rates become too close to each other. 

The essential difference between profile 5, eqs. (9), (10), and the previous one is the 
shape of the decrease along the convective zone. This is enough to ensure a totally different 
quality of the fit. In fact we observe only one very stable minimum in the whole B range 
with an extremely low x 2 (fig- 9): x 2 m i n = 0.054 for Am^ ~ 2.1 x 10~ 8 eV 2 , B ~ 1.45 x 10 5 
G. The stability of this minimum is in sharp contrast with previous cases, as a change of 
15 kG in the peak field implies a change in x 2 by a factor of 2 at the most. An obvious 
conclusion from the comparison with the previous case is that the data prefers a profile with 
an upward facing concavity relative to a downward one along the convective zone. This is to 
be expected, since the low energy sector of the 8 B neutrinos necessarily undergoes a strong 
suppression as the corresponding energy is in the intermediate range, where suppression is 
maximal. Hence in order to ensure a moderate suppression for the whole 8 B flux, its high 
energy part must be suppressed to a much lesser extent implying a strong decrease of the 
field as from the start of the convective zone. It should be noted that such profiles are also 
expected to provide better fits of the high energy part of the recoil electron spectrum in 
SuperKamiokande, as they lead to a weaker suppression of the high energy part of the solar 
neutrino spectrum (see next subsection). 

Profile 6, eqs. (11)-(13), also shows a remarkably good and stable fit. Its general shape 
is roughly the same as the previous one, the only difference being a finite upward slope at 
the bottom of the convective zone. The best fit {x 2 m in = 0.047) is for Am^ ~ 1.60 x 10~ 8 
eV 2 , B ~ 9.6 x 10 4 G (see fig. 10). It is also seen that a change of 11 kG implies a change 
by a factor of 5 in x 2 ■ A second, poorer fit, also shows up in this case: X 2 m %n = 0-65 for 
Am 2 21 ~ 1.3 x 10~ 8 eV 2 , B ~ 2.63 x 10 5 G. 

Finally we examined profile 7. There are four local minima of x 2 m this case, all 
providing fits that are poor and relatively unstable in comparison with the previous ones. We 
have (see fig. 11) : x 2 min = 0-84 for Am 2 21 ~ 1.1 x 10" 8 eV 2 , B ~ 1.26 x 10 5 G, x 2 m in = 0-53 
for Am 2 21 ~ 1.2 x lO" 8 eV 2 , B ~ 1.66 x 10 5 G, x 2 min = 0-71 for Am 2 2l ~ 1.4 x 10~ 8 eV 2 , 
B ~ 2.3 x 10 5 G, x 2 m in = 0-875 for Am 2 2l ~ 1.2 x 10" 8 eV 2 , B ~ 2.7 x 10 5 G. A change 
of 3 kG in the central value of B implies an increase of x 2 by a factor of at least 6. This 
profile is therefore manifestly disfavoured relative to profile 6. It is worth noting that the 
difference between the two lies only in the shape of the downward slope along the convective 
zone. 

3.3 Fits of the recoil spectra and their quality 

In this subsection we calculate the spectrum of recoil electrons in SuperKamiokande for 
each of the magnetic field profiles 1-7. We take in each case the values of (Am 21 , B Q ) that 
correspond to the best fit of the rates as calculated in the previous subsection and evaluate 
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the corresponding x 2 °f the spectrum. We also determine for the cases whose fits and 
stability are the best (profiles 5 and 6) the minimum x 2 an d corresponding (Am| 1; B ). 
The starting point is the ratio of the SuperKamiokande event rate to the one predicted by 
the SSM evaluated for the 18 energy bins (J = 1, 18) 



R = dE e g° dE' e f(E e , g) Jg^ rfE /.„(£?) (P{E)^ + [1 - P(£?)]S|) 
J^dEJ^dEif(E e ,EQJZ UEL) dEfs B (E)^ 



dTdE 



where, f(E e , E' e ) is the energy resolution function which can be found in |j39fl , E e = T + m e , 
and the rest of the notation was defined in section 2. In fig. 12 we plot Rj corresponding 
to profiles 5 and 6 superimposed on data set II. 

As far as the x 2 analysis is concerned there is only one flux now and the values to be 



fitted are the 18 data points given by the SuperKamiokande collaboration (set I ^0] and 



set II [^T| with 708 days and 825 days respectively of data taking). Hence the definition of 



X 2 is now 

18 

X 2 = E (% - R7) (R j2 - R-T) • (35) 

jl,32 = l 

Here the quantities Rj 1 ,Rj 2 denote the vector elements given by (34). Since only one 
neutrino flux is considered here, namely 8 B, in calculating a 2 (ap) the matrix within the 
second summation in equation (32) has now been reduced to its column vector corresponding 
to %\ = %2 = 4 = 8 B. We thus have 

9 

a 2 lj2 (ap) = R h R h £ a 2 4k (A \nX k ) 2 (36) 
fc=i 



where j\, j% = 1, 18. The matrix elements of a(exp) (28) are directly read from PU[ and 



4T| , the matrix o(cs) can be taken to be zero since the uncertainties in the cross sections 
of ve scattering are negligible. 

Evaluating \ 2 f° r the spectrum using the values of (Am^Bo) that correspond to the 
best fit of the rates, we find the results shown in table III for each field profile and data 
set. It is worth noting that the profiles that generate the best and most stable fits (5 
and 6) for the rates are also those that lead to the lowest \ 2 f° r the spectrum. We have 
calculated for these two the minimum x 2 (f° r data set II) and found for case 5 a minimum 
at (3.2 x 10" 8 eV 2 ,1.14 x 10 5 G) (x 2 Pmin = 23.235) to be compared with (see section 3.2) 
(2.1 x 10~ 8 eV 2 , 1.45 x 10 5 G) for the rates. For case 6 (data set II) we found a minimum 
at (2.4 x 10~ 8 eV 2 ,1.12 x 10 5 G) (x 2 Pmin = 23.236) to be compared with (see section 3.2) 
(1.6 x 10~ 8 eV 2 , 9.6 x 10 4 G) for the rates. The stability of these fits should be stressed, as in 
fact the change in the value of x 2 between the two minima (for the rates and spectrum) in 



each of the two profiles only affects the third digit of x 2 ■ F° r this reason we do not include 

(■sp 



the x\ v plots in the paper. 
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We note finally that although the fits for data set II are slightly worse than those for 
data set I, the differences in x 2 are only marginal. Furthermore, one must keep in mind that 
the large uncertainties still present in the high energy (E e > 13 MeV) part of the spectrum 
and the differences in the central data points for sets I and II make it difficult to draw 
any definitive conclusions about the particle-physics solution of the solar neutrino problem 
favoured by the recoil electron spectrum. 



Profile 




xU n ) 


1 


22.5 


24.8 


2 


21.9 


23.9 


3 


22.5 


24.8 


4 


25.9 


29.5 


5 


21.6 


23.5 


6 


21.7 


23.6 


7 


22.3 


24.5 



Table III - The values of x 2 (16 d.o.f.) for the recoil electron spectrum in SuperKamiokande 
evaluated for Atrial, Bo which correspond to the best fits of the rates. Data sets (I) and (II) 
are given in refs. j^j and profiles 1-7 are shown in figs. 1 and 2. Cases 5 and 6, 

besides providing the lowest x% ' s , a ^ so 3™e the lowest x 2 's for the rates and, interestingly 
enough, the most stable fits of all. 



4 Summary and Discussion 



We have examined the status of the RSFP solution to the solar neutrino problem in the light 
of the most recent experimental data on total rates and energy spectrum and the standard 
solar model of BP98 ||. There is independent, although scarce, input information from 
solar physics on the solar magnetic field |!0| suggesting a field that is largest around the 
bottom of the convective zone with a peak up to 3 x 10 5 G in that region. So at the start of 
any RSFP analysis one cannot consider the field profile as entirely free. On the other hand 
previous studies indicate that the common suppression of intermediate energy neutrinos 
and the almost total survival of low energy ones must lead, if neutrinos have a magnetic 
moment, to a profile rising sharply over a relatively short distance along the solar radius, 
decreasing then smoothly towards the surface. Using this input information we investigated 
a general class of profiles with the above characteristics and calculated for each of them the 
values of Am^, B which provide the best fits for the total rates. 

The shape of the survival probability as a function of neutrino energy is obviously 
determinant for the quality of these fits. One important general feature is that the limit of 
this probability as the energy increases is close to 1/2 in the case of RSFP, whereas in the 
case of the small mixing angle MSW solution, which gives the best fit for the rates among the 
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neutrino oscillations solutions, it is close to unity. Hence the important contribution of the 
8 B neutrinos, which survive up to a factor of almost 1/2 and whose energy spread is large, 
is naturally much better fitted in RSFP than in MSW. Thus the fact that the values of Xmin 
for the rates best fits range from 0.036 to 0.59 for 1 d.o.f. (profiles 3 and 4 respectively), 
to be compared with 1.7 (SMA solution) and 4.3 (LMA and vacuum oscillation solutions) 
is not surprising 0. 



The investigation performed for the rates fits has taught us that, besides a fast rise across 
the bottom of the convective zone, the magnitude of the field should decrease faster past 
this bottom and then slower on approaching the solar surface. In other words, the profile 
must show an upward facing concavity. To this end it suffices to compare the fits of profiles 
4 with 5 on one hand and the fits of profiles 6 with 7 on the other. The very low x 2 obtained 
for profile 3 (0.036), with the reverse concavity, is extremely unstable as compared to the 
"best" cases 5 and 6. This fact may be easily understood if one notes that the low energy 
sector of the 8 B neutrinos has to be highly depleted because their energies are close to those 
of 7 Be, CNO and pep neutrinos which must be strongly suppressed by RSFP in order to fit 
the data. On the other hand, the overall suppression of 8 B neutrinos is just moderate. So 
the solar field must ensure a small suppression of their high energy sector. This is achieved 
if the field becomes much weaker soon after the intermediate energy resonance densities. 
Hence an upward facing concavity of the field profile along the convective zone appears to 
be favoured. 

We have also calculated the recoil electron energy spectra in SuperKamiokande and 
performed the corresponding fits for the same class of field profiles. We found for all of them 
except one (profile 4) the correct sign slope against the recoil electron energy. Regarding the 
high energy bins, the central data points show the excess of the observed number of events 
compared with the theoretical predictions, but it should be stressed that the uncertainties 
are much larger in this sector. Hence the general quality of the fits for the spectra is also 
encouraging. We found, with 16 d. o. f., (x%)min ranging from 21.6 to 25.9 (data set I 
40], profiles 5 and 4 respectively) and 23.4 to 29.5 (data set II [H|, also profiles 5 and 4 



respectively). All minima are quite shallow. These values are to be compared with 25.0 
411 ( 24 - x 1!) for the SMA solution, 23.5 g(| for the LMA solution and 17.4 gl| for 



the vacuum oscillation solution []. Future data on the SuperKamiokande energy spectrum 
are thus eagerly awaited, especially regarding the improvement in accuracy of its upper and 
lower ranges. In fact, determining whether the spectrum really rises with energy for E e > 13 
MeV or else the apparent rise is just a statistical fluctuation or an instrumental effect will 
be essential to discriminate between various solutions of the solar neutrino problem. As it 
stands, no theoretical model (oscillations or RSFP) is able to provide a very good fit. 

5 It should be noted that smaller values of Xmin f° r the MSW solutions were obtained in ^o|: Xmin = 0-44 
for the SMA and Xmin = 2 - 7 for tne LMA solution. 

6 The fact that the recoil electron spectrum (including the observed excess of the high energy events) can 
be described reasonably well by vacuum oscillations was first noted in fl4|| and the corresponding solutions 
were found in 122] . 
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It should be noted that reasonably good fits of the spectrum can be achieved in all 
scenarios if one treats the poorly known flux of hep neutrinos, which constitutes the highest 
energy component of the solar neutrino spectrum, as a free parameter j|2], [19|, |2FJ, [|l 



We did not perform such an analysis because the latest SuperKamiokande data shows the 



decreased excess of high energy events |41| and it may therefore be advisable to wait until 
more accurate data has been obtained. 

We have used throughout a value of 10~ n yU,B for the neutrino magnetic moment. Since 
the order parameter is the product \i v B (eq. (17)) and the peak field may be as large as 
3 x 10 5 G, a factor of 2-3 larger than its best fits, the appropriate order of magnitude for 
H v is in fact (3 — 5) x 10~ 12 /xs, a value consistent with most astrophysical bounds |3T|| . 

To conclude, the prospects of RSFP as a possible solution to the solar neutrino problem 
are at present quite encouraging: in fact the fits to the total rates are excellent, far better 



than in any other scenario |18|, ^TJ. As far as the fits for the recoil electron energy 
spectrum are concerned, they are in general slightly better than the small mixing angle 
MSW ones, similar to the large mixing angle MSW ones and worse than vacuum oscillation 
ones, depending of course on the appropriate choice of the field profile. The improvement 
in the data on the spectrum, will be essential in helping to ascertain the most realistic 
theoretical model. 
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Figure captions 



Fig. 1. Solar magnetic fields as functions of the solar coordinate x = r/Rs- Field strength 
is in Gauss. Solid line is profile 1 [eqs. (1) - (3) with xr = 0.7, xq = 0.85], dotted line is 
profile 2 [eqs. (1) - (3) with xr = 0.65, x c = 0.8], dashed line is profile 3 [eqs. (4) - (6)], 
long-dashed line is profile 4 [eqs. (7), (8)]. For x > 0.75 profiles 3 and 4 coincide. 

Fig. 2. Profile 5 [solid line, eqs. (9), (10)] profile 6 [dotted line, eqs. (11) -(13)], profile 7 
[dashed line, eqs. (14) - (16)] . For 0.71 < x < 0.75 profiles 6 and 7 coincide. 

Fig. 3. Neutrino survival probabilities for profiles 1 - 4 as functions of neutrino energy for 
AiUjj and B$ given by their values at the best fit of the rates for each profile. Solid line 
corresponds to profile 1, dotted line, to profile 2, dashed line, to profile 3, long dashed line, 
to profile 4. 

Fig. 4. Same as fig. 3 for the neutrino survival probabilities for profiles 5-7. Solid line 
corresponds to profile 5, dotted line to profile 6 and dashed line to profile 7. 

Fig. 5. Rate fits: values of x 2 1 d.o.f. as functions of the peak value of the field for profile 1 for 
several values of Am^. Left: dashed line corresponds to 7 x 10~ 9 eV 2 , full to 8 x 10~ 9 eV 2 , 
dotted to 9 x 10~ 9 eV 2 . Right: dotted line corresponds to 7 x 10~ 9 eV 2 , full to 8 x 10~ 9 eV 2 , 
dashed to 9 x 10~ 9 eV 2 . Best fit: Am 2 21 = 8 x 10^ 9 eV 2 , B = 1.68 x 10 5 G, xLJd.o.f. = 
0.085. 

Fig. 6. Rate fits: x 2 /d.o.f. as a function of the peak field value for several values of Am^ 
for profile 2. Top left: dashed 1.15 x 10" 8 eV 2 , full 1.2 x 10" 8 eV 2 , dotted 1.25 x 10" 8 eV 2 . 
Top right: full 1.3 x 10" 8 eV 2 , dotted 1.35 x 10" 8 eV 2 , dashed 1.40 x 10" 8 eV 2 . Bottom left: 
dotted 1.15 x 10" 8 eV 2 , full 1.20 x 10" 8 eV 2 , dashed 1.25 x 10~ 8 eV 2 . Bottom right: dotted 
1.15 x 10~ 8 eV 2 , dashed 1.20 x 10~ 8 eV 2 , full 1.25 x 10" 8 eV 2 . Best fit: Am 2 21 = 1.20 x 10~ 8 
eV 2 , B = 1.23 x 10 5 , xLJ d.o.f. = 0.100. 

Fig. 7. Rate fits: x 2 / d.o.f. as a function of the peak field (profile 3) for several values of 
Am^. Left: dotted 1.1 x 10" 8 eV 2 , full 1.2 x 10" 8 eV 2 , dashed 1.3 x 10" 8 eV 2 . Right: 
dot-dashed 1.1 x 10~ 8 eV 2 , full 1.2 x 10~ 8 eV 2 , dotted 1.3 x 10~ 8 eV 2 , dashed 1.4 x 10~ 8 
eV 2 . Best fit: Am 2 21 = 1.2 x 10~ 8 eV 2 , B = 9.54 x 10 4 G, xliJd.o.f. = 0.036. 

Fig. 8. Rate fits: x 2 /d.o.f . as a function of the peak field (profile 4) for several values of 
Am 2 !. Left: dashed 1.1 x 10" 8 eV 2 , full 1.2 x 10~ 8 eV 2 , dotted 1.3 x lO" 8 eV 2 . Right: dashed 
1.5 x 10- 8 eV 2 , full 1.6 x 10~ 8 eV 2 , dotted 1.7 x 10" 8 eV 2 . Best fit: Am 2 x = 1.6 x lO" 8 eV 2 , 
B = 1.7 x 10 5 G, xLJd.o.f. = 0.59. 

Fig. 9. Rate fits: x 2 /d.o.f. as a function of the peak field for profile 5 for several values of 
Am 2 x . Dotted line 1.9 x 10~ 8 eV 2 , dashed 2.0 x 10~ 8 eV 2 , full 2.1 x 10~ 8 eV 2 , dot-dashed 
2.2 x 10" 8 eV 2 . Best fit: Am 2 2l = 2.1 x 10" 8 eV 2 , B = 1.45 x 10 5 G, X 2 m iJ d -°-f- = 0.0547. 

Fig. 10. Rate fits: x 2 1 'd.o.f. as a function of the peak field for several values of Am 21 
(profile 6). Left: dotted 1.55 x 10" 8 eV 2 , full 1.6 x 10" 8 eV 2 , short dashed 1.65 x 10" 8 eV 2 , 
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long dashed 1.70 x 1(T 8 eV 2 , dot-dashed 1.75 x 1(T 8 eV 2 . Right: dotted 1.2 x 1(T 8 eV 2 , 
full 1.3 x 1(T 8 eV 2 , dashed 1.4 x 1(T 8 eV 2 . Best fit: Am 2 21 = 1.6 x 1(T 8 eV 2 , B = 9.6 x 10 4 
G, xLJd.o.f. = 0.0477. 

Fig. 11. Rate fits: x 2 /d-o.f. as a function of the peak field (profile 7) for several values 
of Am 2 21 . Top left: dotted 1.0 x 10~ 8 eV 2 , full 1.1 x 10~ 8 eV 2 , dashed 1.2 x 10~ 8 eV 2 , 
dot-dashed 1.3 x 10~ 8 eV 2 . Top right: dotted 1.1 x 10~ 8 eV 2 , full 1.2 x 10~ 8 eV 2 , dashed 
1.3 x 10~ 8 eV 2 . Bottom left: dotted 1.0 x 10" 8 eV 2 , short dashed 1.1 x 10^ 8 eV 2 , long 
dashed 1.2 x 10" 8 eV 2 , dot-dashed 1.3 x 10~ 8 eV 2 , full 1.4 x 10" 8 eV 2 . Bottom right: dotted 
1.0 x 10" 8 eV 2 , dot-dashed 1.1 x 10" 8 eV 2 , full 1.2 x 10~ 8 eV 2 , short dashed 1.3 x 10~ 8 eV 2 . 
Best fit: Am 2 x = 1.2 x 10" 8 eV 2 , B = 1.66 x 10 5 G, xLn/d-o.f. = 0.53. 

Fig. 12. Spectrum fits: for profiles 5 (top) and 6 (bottom) we show the predicted electron 
recoil spectrum as a function of total electron energy superimposed on data set (II) (825 
days), ref. ^T|, with Am^ and B given by their values corresponding to the rates best fits 
(see figs. 9 and 10). For profile 5 x%/d.o.f. = 23.5/16, and for profile 6 x%/d.o.f. = 23.6/16. 
For both x% /d.o.f. = 23.2/16 with values of Am| 1; So that are slightly different from 
the best fit ones for the rates (see the main text). 
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